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^ Abstract 

A: 

We present the results for three-loop gauge field anomalous dimensions in the 
SM calculated in the background field gauge within the unbroken phase of the 
model. The results are valid for the general background field gauge param- 
eterized by three independent parameters. Both quantum and background 
fields are considered. The former are used to find three-loop anomalous di- 
mensions for the gauge-fixing parameters, and the latter allow one to obtain 
the three-loop SM gauge beta-functions. Independence of beta-functions of 
gauge-fixing parameters serves as a validity check of our final results. 



1 Introduction 



In spite of the fact that the Standard Model has many unsatisfactory aspects 
Nature still does not allow us to find some solid evidence for the existence of 
a more fundamental theory with new particles and/or interactions. Due to 
the joint efforts of both experimentalists and theoreticians we are about to 
enter the only unexplored part of the SM and unveil the mechanism of elec- 
troweak symmetry breaking. According to the recent experimental results, 
there is strong evidence for the existence of the Higgs boson, the last missing 
ingredient of the SM spectrum |lll2]- 

The mass of the higgs seems to be located at the boundary of the so- 
called stability and instability regions in the SM phase diagram. This fact 
implies that the SM can be potentially valid up to a very high scale (e.g.. 
Plank scale). 

In this situation, it is important to know how the running SM parame- 
ters evolve with energy scale. The analysis of high energy behavior is usu- 
ally divided into two parts. The first one is the determination of running 
MS-parameters from some (pseudo)observables. This procedure is usually 
referred to as "matching". The second one utilizes renormalization group 
equations (RGEs) to find the corresponding values at some "New Physics" 
scale. In order to carry out such an analysis consistently one usually use 
(L — l)-loop matching to find boundary conditions for L-loop RGEs (see, 
e.g., [3j). It is worth pointing that the advantage of the minimal-subtraction 
prescription lies in the fact that one needs to know only the ultraviolet (UV) 
divergent part of all the required diagrams. The latter has a simple polyno- 
mial structure in mass and momenta (once subdivergences are subtracted). 
Due to this, MS beta functions and anomalous dimensions can be relatively 
easily extracted from Green functions by solving a single scale problem with 
the help of the so-called infrared rearrangements (IRR) 

One- and two-loop results for SM beta functions have been known for 
quite a long time |5l El III El O [TOl HI [El US E] and are summarized in [15] . 
Until recently, three-loop corrections were known only partially [T71 [T5| 
da EDI ED- 

Having a well tested method for calculation of three-loop renormalization 
constants [22l [231 121] and an experience in the calculations in the Standard 
Model and its minimal supersymmetric extension [251 EEl |27] we are planning 
to perform the calculation of all renormalization group coefficients in the 
third order of perturbation theory extending the results of Refs. [IHl EHl [29] 
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to one more loop. 

In this paper, we present our first step in this direction: the results for 
three-loop anomalous dimensions of the SM gauge fields. Since we are only 
interested in UV-divergences for the fields and dimensionless parameters, we 
do not consider the effects related to spontaneous breaking of electroweak 
symmetry and, as a consequence, can neglect all dimensionful parameters 
of the model. Moreover, we made use of the background-field gauge (BFG) 
(see, e.g., Ref. [30] and reference therein) to carry out our calculation. In 
this gauge, due to the simple QED-like Ward identities involving background 
fields, one can easily obtain expressions for the beta-functions by considering 
the two-point functions with external background particles. 

During the work on this project a few papers on the same topic appeared 
[311 [32] (gauge couplings) and [33] (Top Yukawa and higgs self-interactions). 
Since the authors of [32] carried out a similar calculation, let us mention that 
our setup differs from that used in Ref. [31] in several aspects. 

Firstly, for the diagram generation we made use of FeynArts [31], as it 
includes well-tested model files for the SM. Since the diagrams are evaluated 
with the help of the MINCER package [3B], a mapping to the MINCER notation 
for momenta is required. This problem was solved by hand with the help 
of the DIANA [36] topology files which were prepared during our previous 
calculations p2] . Based on these files a simple script was written which allows 
one to perform the mapping between the FeynArts and MINCER notation. 

Secondly, we not only consider Landau BFG in the broken phase of the 
but also choose to work within the unbroken SM in a general BFG. The 
full dependence of the diagrams on the electroweak gauge-fixing parameters is 
retained and the corresponding renormalization is taken into account. And 
lastly, since the unbroken SM in BFG is not implemented as a FeynArts 
model file, we are forced to use a package like FeynRules [57] or LanHEP [55] . 
Due to the fact that the authors are more accustomed to the latter, it was 
chosen to generate the required Feynman rules from the Lagrangian. 

The paper is organized as follows. In Section [2] we introduce our nota- 
tion and present a brief description of the unbroken SM quantized in the 
background-field gauge. Section [5] describes the details of our calculation 
strategy. Finally, the results and conclusions can be found in Section HI Ap- 
pendix contains all the expressions for the considered renormalization con- 
stants. 

^ As a cross-check of our main calculation within the unbroken SM. 
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2 The Standard Model in the unbroken phase. 
The background-field gauge. 

Let us briefly review the Lagrangian of tlie SM in tlie background-field gauge. 
We closely follow [32] albeit the fact that we introduce background fields only 
for gauge bosons. Moreover, as it was mentioned in Introduction, we neglect 
all the dimensionful couplings (i.e., mass parameters). 
In our calculation we use the Lagrangian of the form 

C = Cg + U + Cf + jCgf + Cfp. (1) 
Here £g is the Yang-Mills part 



1 1 . . 1 

-G" G" - -W W - - 
4 4 4" 

= d^Gl — d^G^j^ + gsf"''"'G^^Gl, (3) 
K'^ = d.Wl, - duW; + g2e'''Wj^Wt. (4) 
= d^B, - d,B^, (5) 

where G» = + G^l (a = 1,...,8), W^, = W; + W^, (z = 1,2,3), and 
B^ = B^ + B^ are gauge fields for SU(3), SU(2) and U(l) groups. By 
V = {G, W, B) we denote quantum fields, and V = {G, W, B) is used for 
their background counterpart. The corresponding gauge couplings are Qs, g2, 
and Qi. The group structure constants enter into the commutation relations 

^ya^yfcj _ .^abcyc^ [r\ r^] = ze*^' V, (6) 

with T" = A"/2 and r* = (7^/2 being color and weak isospin generators. 

The covariant derivative acting on a field which is charged under all the 
gauge groups looks like 

D, = d,- gsT^Gl - g^rW; + g,^B,. (7) 

If a field is not charged under either group, the corresponding term is omitted. 
With the help of the covariant derivative one can write the following Higgs 
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and fermionic parts of the Lagrangian: 

Cn = {D,<^)HD,<^)-^{^^^)\ (8) 

i=l,2,3 



- E 

ij=l,2,3 

where indices i,j = 1,2,3 count different fermion famihes, A and Yu^d,i are 
the higgs quartic and Yukawa matrice^, respectively. The left-handed quarks 
Qg = (""35 (^g)^ leptons Lg = {ug, Ig)^ form the SU(2) doublets while the 
right-handed quarks (u^, d^) and charged leptons are the singlets with 
respect to SU(2). The Higgs doublet $ with Yvv = 1 has the following 
decomposition in terms of the component fields: 



.'^.^^^[^-p'y (10) 

Here a charge-conjugated Higgs doublet is introduced with Yw = —1. 
The gauge-fixing terms are introduced only for quantum fields 

1 1 1 

^Kg ^Kw ^Kb 



with 



G\v = d,W; + g2e'^'Wj,Wl 
GB = d,B^. (12) 



rk 
■ At ' 



The ordinary derivatives are replaced by covariant ones containing the back- 
ground fields. Due to this, the invariance of the effective action under back- 
ground gauge transformations is not touched by introduction of (ITTl) . 
The Fadeev-Popov part of the Lagrangian is given by 

/^FP = (13) 



^In the actual calculation the diagonal Yukawa matrices were used. However, the result 
can be generalized with the help of additional tricks (see Secl3]and Ref. 32 ). 
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where = {G,W,B), and 5Ga/S6^ is the variation of gauge-fixing func- 
tions fll2l) under the following infinitesimal quantum gauge transformations 

^B^ = O^Ob . (14) 

It should be stressed that covariant derivatives in f|T4l) involve the sum of 
quantum and background gauge fields V = V + V. The corresponding back- 
ground transformations are obtained from fll4p by the replacement V V. 

The Feynman rules for the model described by the Lagrangian ([1]) were 
generated with the help of LanHEp[§ [38j . 

It is worth mentioning here that our problem does not require the intro- 
duction of U(l) ghosts cb,cb and background 13 fields. This is due to the 
fact that the latter has the same interactions as its quantum counterpart B 
and the former decouples from other particles. Nevertheless, we keep them 
in our LanHEP model file to allow for possible generalizations to non-linear 
gauge-fixing as in Ref. [39] . 



3 Details of calculations 

Due to the gauge invariance of the effective action for the background fields, 
QED-like Ward identities can be derived. The latter can be used to prove 
the following simple relations: 

Z,, = Z7^/^ z = l,2,3 (15) 

with Zy^ and Zg. being renormalization constants for background fields = 

{B'^,W'^,G^) and SM gauge couplings gi = {gi,g2,gs), respectively. 

Since we keep the full dependence on the gauge-fixing parameters dur- 
ing the whole calculation, we also need to know how = {^b,^w,^g) sjce 
renormalized. Again, due to the Ward identities, the longitudinal part of the 
quantum gauge field propagators does not receive any loop corrections. As 
a consequence, the following identities hold: 



■^LanHEP 3.1.5, which was used by the authors, produces a wrong sign for the combi- 
nation f^bcjdec ^■^^Y[lJg export to the FeynArts model files. A new version with a fix is 
scheduled for November 2012. 
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Here Z^. stands for the renormalization constants for the gauge-fixing pa- 
rameters. The quantum gauge fields Vi are renormahzed in the MS-scheme 
with the help of Zy_ . It is clear from f lTSj) and f|T6l) that to carry out the cal- 
culation, one needs to consider gauge boson self-energies for both quantum 
V and background V fields. 

For calculation of the renormalization constants, following [12] (see also 
[HI in Sn]), we use the multiplicative renormalizability of the corresponding 
Green functions. The renormalization constants Zy relate the dimension- 
ally regularized one-particle-irreducible two-point functions Fy^Bare with the 
renormalized one Fy^en as: 

rv,Ren ( ^^'^i ] = hmZy ( -, Oj J FyBare (Q^, C-i.Bare; ; (17) 

where Cj^are are the bare parameters of the model. For convenience, we 
introduce the following notation, which is closely related to that used in 
Ref. 



a. = (-^, ^G, iw. ic 

so we treat the gauge-fixing parameters along the same lines as couplings. 
Moreover, in the renormalization group analysis of the SM one usually em- 
ploys the SU(5) normalization of the U(l) gauge coupling which leads to an 
additional factor 5/3 in f lTSj) . 

The bare parameters are related to the renormalized ones in the MS- 
scheme by the following formula: 



flfc.Bare/^ 



oo ^ 

-'''^' = Z^M^^) = au + Y.ct^-. (19) 



n=l 



where pk = I for the gauge and Yukawa constants, pk = 2 for the scalar 
quartic coupling constant, and pk = for the gauge fixing parameters. In 
order to extract a three-loop contribution to Zy from the corresponding self- 
energies, it is sufficient to know the two-loop renormalization constants for 
the gauge couplings and the one-loop results for the Yukawa couplings. This 
is due to the fact that the Yukawa vertices appear for the first time only in 
the two-loop self-energies and the higgs self-coupling enters into the result 
only at the third level of perturbation theory. 
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The four-dimensional beta- functions, denoted by are defined via 

dai{fi,e] 



(20) 



e=0 



(iln/i^ 

Here, again, stands for both the gauge couphngs and the gauge-fixing. 

Given the fact that the bare parameters do not depend on the renormal- 
ization scale the expressions for /3j can be obtained [13] by differentiation of 
( !T9|) with respect to In/i^: 



OO ^ ■ 



n=l 



dc 



(n) 



n=l i 



(21) 



Taking in account only the leading order of the expansion in e: 

/3k = y^^Piai- 



dc^' (1) 



(22) 



In MS-like schemes the renormalization constants for the Green functions 
may be expanded as 



(23) 



k=l 



Differentiating fl23p with respect to In /x^ we simply get all-order expression 
for anomalous dimensions: 



7r 



-P 



,d\nZr 



dZr 



(24) 



It turns out that the above expression is finite as e — )■ so 



7r 



(25) 



The advantage of (!?!]) and comes from the fact that it provides us 
with additional confirmation of the correctness of the final result since beta 
functions and anomalous dimensions extracted directly from fl2T]) and 
are finite for e ^ only if c^"'' satisfy the so-called pole equations [H], e.g.. 



d 



(26) 



In order to calculate bare the two-point functions for the quantum and 
background fields, we generate the corresponding diagrams with the help of 
the FeynArts package J3^. It is worth pointing that we use the Classes level 
of diagram generation which allows us to significantly reduce the number of 
generated diagrams since we do not distinguish fermion generations. The 
complexity of the problem can be deduced from Table [T] that shows how the 
number of the FeynArts generated diagrams increases with the loop level. 
Clearly, the presented numbers are an order of magnitude less than those 
given in Table I of Ref. [32], which somehow demonstrate the advantage of 
our approach. 

The number of the SM fermion generations is introduced by hand via 
counting fermion traces present in the generated expression for a diagram 
and multiplying it by no- We separately count fermion traces involving the 
Yukawa interaction vertices and multiply them not by no but by ny. This 
allows us to use the following substitution rules (c.f., [32]) to generalize the 
obtained expression to the case of the general Yukawa matrices 

nY[au,ad,ai] [3^«, 3^d, 3^/] , 
nY[al,al,a^] [yuu,ydd,yii], 



n 



Y 



Kaia^] ^ [ylylyf]-^ 



n 



Y 



[yuyd,ydyi,yuyi], 

nyauad (27) 



where 



and 



^""-^^^ ^^-^6^^' ^'-^6^' ^^^^ 



ti Y^Y^Y^Y^ ^, _trYdY}YdYj 



~ (167r2)2 ' ^'"'^ (167r2)2 ' 

^^g^2^2 ' ^« (16^2)2 • ^^^^ 



A comment is in order about the last substitution in fl271) . It turns out 
that yud is the only combination of up- and down-type Yukawa matrices, 
which can appear in the result for the three-loop gauge-boson self-energy 
within the SM. This can be traced to the following facts: 1) in the unbro- 
ken SM all the particles are massless so that chirality is conserved; 2) the 
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Broken 


1 


2 


3 


Unbroken 


1 


2 


3 




10 


339 


21942 




11 


389 


36647 


z 


9 


281 


19041 




11 


371 


36103 


A 


7 


218 


14426 


B,B 


6 


214 


20144 


ZA 


7 


236 


16120 


G 


4 


73 


4183 


G 


4 


67 


3287 


G 


4 


66 


4060 


Total 


37 


1141 


74816 


Total 


36 


1113 


101137 



Table 1: Number of self-energy diagrams with external gauge fields, gener- 
ated by FeynArts in the broken and unbroken SM, at one, two, and tree 
loops. 



Yukawa interactions flip the chirality of the incoming fermions; 3) there is 
no right-handed flavour changing current coupled to a SM gauge field. As a 
consequence, combinations like 

l^yytyyt irYYW,,Y'^ 

tiYuy.Yuya (30) 



(167r2)2 (IGtt 

which require at least two chirality-conserving transitions between right- 
handed up- and down-type quarks, do not show up in the result. 

This type of counting is performed at the generation stage. A simple 
script converts the output of FeynArts to DIANA-like [3S] notation and iden- 
tifies MINCER topologies. This allows us to use the FORM |l2] package COLOR 
[43j to do the color algebra and MINCER [35] to obtain the e-expansion of di- 
agrams. It is worth pointing that the expressions for all SM gauge couplings 
exhibit explicit dependence on number of colors A^^^ which stems from the 
fact that we have to sum over color when there is a (sub) loop with external 
color singlets coupled to quarks. 

Having in mind the cancelation of gauge anomalies within the SM we use 
a naive anticommuting prescription for dealing with 75, so that all the Dirac 
traces involving one 75 vanish. 

4 Results and conclusions 

Here we present the results of our calculations in the form of the SM gauge 
beta-functions and anomalous dimension of the gauge-fixing parameters. 
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From fllSI) and (fTB]) it is clear that anomalous dimensions of the background 
fields are connected with the corresponding gauge coupling beta-functions 



7b = -/5i/«1' 7w = -/32/a2, 7g = (31) 
and for the quantum fields we have 

7B = /3«s/Ci3, lw = P^wl^w, lG = PiGliG- (32) 

The corresponding renormalization constants can be found in the Appendix. 

At the end of the day, we have the following expressions for the beta- 
functions: 

fllN^ 3\ 1 



,2 



/137aiiV, 81ai a2 9a2 lla.C^ A^'^ 



+ «i 1 I — ^7^7^ — + + + — + 



900 100 20 20 15 



9ai 9a2 N^y^ llN^y^ 3^^ 
50 10 6 30 2 



ah no 



imial 981 Nc 27aia2 

18000 2000 1200 400 



137 ^ al N, 27al 1 ^ 1463 ^ 

900 ^ ' ^ 45 ^ 10 20 540 " 

^^a^C'^N\ r? /_16577a2A^2 23870^ A^^ _ 8910^ WalNl 

30^" ^ + 486000 9000 2000 720 

lla^A-^ 11 a2 242 2^ \ 4890^ 783ai 02 27aiA 

:a CpTp Ac + H — — h 



72 80 135 ' 7 8000 800 50 

1267ai A'cJ^d 2827aiA'c 2529ai3^i 3401 al 9a2A 

2400 2400 800 ^ 320 ^ 10 

437a2Ar,3;, 157a2Ar, 1629a23^i 17 ^ 29 ^ 

a=CpAcJ/w cIsCp AcJ/u 

160 32 160 20 ^ ^ ^^"^ 20 

_ 9A2 17Ar^^ ^ V ^ ^c^' ^ 157Ar,y, y, 61Ar,y,, 

5 ^ 120 + 60 ^^'^"^ 120 ^ 60 ^ 80 

^ 199A^c yiyu ^ N,y^d ^ 113 N,y^^ 993f 261_3^\ 

60 + 8 + 80 + 40 + 80 r ^"^"^^ 
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o 43 
/32 = a,(nG (- + -)-- 



+ «9 'T'G K a5<-^ F^c 

V 60 20 12 12 ^ ^ 

^3ai 259a2 TVc A^cX 3^^ 



10 6 2 2 2 

,2/ 287a2Ar, 91a2 13aia2 AT, 39aia2 1 „ 

+ "^l^<-^60^ - 400 + ^40— + ^"60^^ ""^^'^^ 
mSal iVc 1603a^ 13 ^ 133 2^ ^ 1 2^2 



+ nl\ - 



121al Nl _ na\Nc _ 33 al _ AlhalNl _ 415 ajA^c _ 415a| 
32400 1800 400 432 216 432~ 



22 2^ ^ , A , 163 al , 561aia2 , 3ai A 533aiA^e3^d 593 aiA^c3^„ 



9 ^ ^ ^ y 1600 160 10 480 480 

_ 51ai yi _ 667111ai 3a2 A _ 243a2A^c3^d _ 243 a2A^c3^„ _ 243a2 yi 
32 1728 ^ 2 32 32 32 

- '[asCFN, ya - \asCFN,y^ - 3 A^ + + ^Nj y^yu + 



16 4 8 16 

16 



+ ^ + ^1. (34) 



^T^nc 11 Ca 



3 3 



+ I "-G I — h 3 a27V H g h 8 asCpTp 



^( (- 



- ^— - 4 T^:;, - 4r^3;„ J + a, [no ^ - — 

^ 22 11 , 241air^ 
+ — aiasCAiF - —ai asCpiF H — 

+ 6a2asCATF - 3a2asCF Tp + ^^^a^C^Tp + gICaCfTf - 4:alCp Tp 

11 



+ nr[ a<)iF^c ^^(-'A-Lp 

^\ 8100 300 12 2 12 27 ' ^ 

176 2^ 89aiT^3^d 101 aiT^X Q^a2Tp 93a2rpX 

-— a.GW^J 20 4 4 

2857 

- —r a^Cl - 2Aa,CATFyd - '24. asC^TFy^ - ^asCpTp y^ - ^asCpTpy^ 

+ itf NX + uTpN^y^yu + 7 TfNX + 7TFydyi + 9 TFy^ + ^TFylyu 

- qtf yud + QTFyuu] , (35) 



/ ai / 3ai llaiiVc\\ 

50 10 6 2 30 

/ 81a? 9aia2 137a?iVc aiaziVc H ^ 
+ nc[ aiasGpA'c 

V 100 20 900 20 15 J J 

/ 489a? 783a2a2 3401aia^ _ 54a?A _ 18aia2A 144aiA2 
^ '^^ V 8000 800 320 25 5 ^ 5 

/981a? 27a?a2 27aia2 1697a?iVc a?a2iVc 1 . , 
+ 71(7 1 1 1 QiaiNf. 

V 2000 400 10 18000 1200 45 ^ 

137 2 ^ AT 1 ^ 1463 2^ ^ AT 11 2^2 AT \ 

900 ^ 20 540 ^ 30 ' 7 

2 /^891a? llaiai ^ 2387a?7V^ 11 2 a, ^ 242 . 

+ GiaiNc H aiaXpTpA/c 

^V 2000 80 9000 72 ^ 135 ^ 

16577a?Ar2 11 .A 1267alNcyd 437 

H ^ + aialN^ + ^ aia2iVe3^d 

486000 720 ^ V 2400 160 

^1"^ ^; /\r2^;2 61ai7V,yrfrf 2529a?y, 16290102^; 
+ ^^a.a.CFN^yd - ^^a,N,ya ^ + + igQ 

157 99aiyf 261aiya , 2827a?iVey„ , 157 
- -^a, ^ + ^^QQ + ^aia2iV.X 

29 59 199 101 

+ -a,asCFN,y^ - -a,N!y,y, - —a,N,yiy^ - —a^Nlyl 

20 60 60 IzO 



8 80 



12 



/25a2 f a2 QiN, 



V 6 ■ ""V 3 3 )) 

/ 3aia2 llSol llal^w 2^2 Z' 3aia2 

1^-^ + ^ 2 + 



20 



4 60 4 
/ 163a?a2 33aia2 143537a| 6aia2A \2 
+ ^^(,-^600 32- + ^76 ^-6a2A + 48a2A 

_ 315ai^w _ 33ai^^ _ 7ai^^ 2 /33a?a2 185ai 77aja2Nc 
8 4 ^ ^ V 400 ^ 144 ^ 1800 

185aiiVe 22 2^ ^ 121a2a2A^c^ 185alN^\ 533 

72 ^ 9 ^ ^ ^ ^ ^ 32400 144 y/ 480 

32^2 '' '^ |a2a5Cj.iVc3^d - ga2A^c3^d 5^ + ^2 

79alyi 5 5a23^f 19a2y» ^ 593 ^ 79 ^ 

+ ^a2a,C^7V,X - ^-a2N'^y,y^ - ^a2iVe3^,X - la2N'^yl - ^^H^ 

19Q2A^cX. 9 , 3aiC(3) ^ 3 . 3 3 . 
Y^aia2C(3) + 6a2CiyC(3) - 2«2WC(3) 

/91a? a2 Qaial 7025a^ ^ 287a?a2iVc 2 

400 5 144 3600 15 ^ 

7025a^iV 1 133 

—7—^ + —aia2asCFNc + SalagCpNc - —a2a1CACFNc 

144 60 36 

+ ^a2alC^FNc + 2aUwNc - ^a,alC{3) + 9a^C(3) - ^aialN,C{3) 

+ 9ai7VeC(3) - 12alasCFN,C{3)^ ^ , (37) 



/13asCA OsCaCg 8asTFraG\ 
^^--^^^~6 2 3~) 



+ Cg - y a^Cl^G - ^a^Clel + ^a^TFya + 4a,r^3^„ 



13 



11 



11 

15 



aia^Tp — 3a2asTp — lOa^C^Tp — Sa^CpTp 



/9965a3c3 167 



288 



32 



''"^ ' 3/~f3 3/^3 /-2 ' 3/^3 /-3 2 



7 

32' 



A21 2 ^ 
V300 ^ 



304 



176 



2 l?,?,lalasTpNc 11 2 



12^ gsAF 9^^^ 8100 



93 

T 



25 

Y 



+ ^aiasTi7 3^d + ^a2asTpyd + -^alCATpyd + GalCpTpyd - lasTpN^yl 



101 93 25 

- 9asTpydd - 7asTpydyi + —aittsTpyu + —a2asTpyu + —alCATpyu 

+ QalCpTpyu - UasTpN.ydyu - lasTpyiy^ - TasTpN.yl + QasTpy^d 
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- Qa^Tpy^^ - i-alCla^) - -XClicm - T^a^C^CcClS) + nci i^afasTp 



3,^3 £-2 



16 



4 



16 



13 



60 



1 ^ 241 2 ^ 319 2^ ^ 87 2r. ^ 911 3^.2^ 

+ -^0'ia20'sTp — —a2asTp + — — aia^C^iTF + —a2agCATp — a^C^Tp 

zu iz izu o y 

+ ^aia^C^Tp + 3a2alCpTp - ^alCACpTp + AalClTp + 4a3CleGTp 
15 9 

22 

- — aia^C^rpC(3) - 18a2a2C^TpC(3) + 3QalClTpC{3) 





- 48a^CACFTFC(3) 



(38) 



With the help of substitutions = A^-^ = 3, = 4/3, = 1/2, 
y^ = trT, yd = tiB, y = tiL, ydd = tr{B^), yuu = tr(f2), yu = tr(f2), 
and yud = tr TB it is possible to prove that the expressions presented above 
coincide with the results for the gauge beta functions obtained in Ref. [3T] . 

As a consequence, one can be sure that the three- loop renormalization 
group equations obtained for the first time in Ref. [31] are correct and con- 
firmed by an independent calculation. It is also worth mentioning that the 
obtained results can be used not only for the analysis of vacuum stability 
constraints within the SM (see, e.g., [HJ HSl |16]) but also, e.g., for very 
precise matching of the SM with its supersymmetric extension since the cor- 
responding three-loop renormalization group functions are already known 
from the literature [TTl HH] . Moreover, the leading two- loop decoupling cor- 
rections for the strongest SM couphngs are also calculated within the MSSM 
in Refs. [251 ill |50l [26] . 
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A Renormalization constants 

Here we present the results for the renormalization constants from which 
the anomalous dimensions and beta-functions were extracted. It should be 
pointed out that the coefficients of the e-expansion satisfy the pole equa- 
tions (126|) . The correspondingexpressions together with the results for beta- 
functions can be found onlinal in the form of Mathematica files. 
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